function [ f ] = fitFunction( x, xdata, sigma )

omega = 4*pi;

epsilon = sqrt(omega*sigma./(xdata*pi));

f = 1/15 + epsilon*(1/3)*1/x;

end

